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ABSTRACT 

In this paper, we present systematic studies on the B-, R- and Kg-hand luminosity- 
metaUicity {L-Z) relations for a set of metal poor, blue compact dwarf galaxies. Metal- 
licity is derived by using both the empirical N2 and the direct Tg methods. Our work 
reconciles contradictory results obtained by different authors and shows that the L-Z 
relationship does also hold for blue compact dwarf galaxies. The empirical N2-based 
slope of the L-Z relation, for each photometric band, is consistent with the Tg-based 
one. We confirm that the slope of the L-Z relation is shallower in the near-infrared than 
that in the optical. Our investigations on the correlations between the Lb-Z relation 
residuals and different galactic parameters show that the star formation activities could 
be a cause of the large scatter in the optical L-Z relationships, whereas the internal- 
absorption might be another possible contributing factor. 

Subject headings: galaxies: dwarf - galaxies: spectroscopy - galaxies: starburst- galax- 
ies: abundances 



1. Introduction 

Since the seminal work of Lequeux et al. (1979), the relationship between luminosity {L) and 
metallicity {Z; hereafter L-Z relation), two crucial parameters for understanding the behavior of 
galaxies, has been extensively studied (e.g. Skillman et al. 1989; Tremonti et al. 2004). Metallicity 
reflects the gas reprocessed by stars and any exchange of gas between environments and galaxies. 
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while stellar mass (i.e. luminosity for constant mass-to-light ratio) reflects the amount of gas 
locked up into stars. Therefore, the relationship between luminosity and metallicity can provide 
an important constraint on the models of galaxy formation and evolution that attempt to account 
for the chemical evolution of the system (see e.g., Prantzos &: Boissier 2000; Boissier et al. 2003; 
Qian & Wasserburg 2004). Additionally, the L-Z relation could be also used to search for the most 
metal-poor galaxies and to derive, at least to a first approximation, the distance of a galaxy from 
its metallicity and vice versa. 

A well-defined correlation between the blue luminosity (Lb) and the metallicity for dwarf 
irregulars (dis), spirals and ellipticals has been established by various authors (e.g., Garnett & 
Shields 1987; Skillman et al. 1989; Zaritsky et al. 1994; Melbourne & Salzer 2002; Lee et al. 2003; 
Lamareille et al. 2004; Tremonti et al. 2004; Salzer et al. 2005; van Zee & Haynes 2006), over ~10 
mag in luminosity and 2 dex in metallicity. However, some studies do not support these results 
for all types of galaxies. Specially, several contradictory results regarding the existence of an L-Z 
relationship of dwarf galaxies have been presented over years. Hidalgo-Gamez & Olofsson (1998) 
presented a weaker Z-Z relationship than previously thought, between Mb and 0/H for dIs. Hunter 
&; Hoffman (1999; hereafter HH99) found that the relationship between Mb and 0/H for Im, Sm 
and blue compact dwarf galaxies (BCDs) has a very large scatter. Another question is whether the 
L-Z relationship for dwarf galaxies (if it is true) exists in a similar manner as for massive galaxies. 
Using a sample of more than 500 star-forming galaxies, Melbourne & Salzer (2002) found that the 
L-Z relation for giant galaxies has less scatter and is steeper than that for dwarf galaxies (also see 
Salzer et al. 2005). Using ~1000 individual spectra of Hll regions in 54 late-type galaxies, however, 
Pilyugin et al. (2004) found that the slope of the L-Z relationship for spirals (—18 < ALb < —22 
mag) is slightly shallower than the one for dIs (— 12 < Mb < — 18 mag), which is consistent with 
the result for blue compact galaxies (BCGs; Shi et al. 2005). 

BCDs are a group of extragalactic objects that are spectroscopically characterized by a faint, 
blue optical continuum accompanied, in most cases, by intense emission lines. BCDs are small 
galaxies with low metallicities (1/50 < Z < 1/3 Zq; HH99) and have dramatically different prop- 
erties compared to normal dwarf galaxies (Zwicky 1966; Gil de Paz et al. 2003). These galaxies 
have bluer colors than ordinary dIs and their surface brightness is much higher. Accumulated ob- 
servational evidence over the recent years provided more details on the unique properties of these 
galaxies (for a review see Kunth & Ostlin 2000). Having a relatively low metallicity, BCDs are at 
an early epoch of their evolution, making them similar to samples of the distant, more massive pro- 
togalaxies, thus allowing us to study the star formation and chemical enrichment in an environment 
likely to be similar to that in the early universe. 

Although it is thought the L-Z relationship should exist in almost all types of galaxies, it is not 
clear yet whether this relationship holds for BCDs alone (HH99, Fig. 7; Hopkins et al. 2002, Fig. 
7). Vaduvescu et al. (2007) show that there exists a good correlation between the metallicity and 
the Kg luminosity for BCDs. However, results from these authors may suffer large uncertainties 
from the limited size of their BCD sample (13, 22 and 14 galaxies in HH99, Hopkins et al. and 
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Vaduvescu et al., respectively), and/or the selection criteria for BCDs, as well as the inconsistence 
in the determination of metallicity for different objects. Therefore, further works are needed to 
disentangle these two apparently contradictory results. 

In this paper, we present our systemic studies on the L-Z relationships for BCDs from the 
optical to the near-infrared (NIR), i.e. at the i?-, R- and -fC^-band. Our BCD sample is derived 
from the data in Gil de Paz et al. (2003), and it contains ~ 70 galaxies (40 of which have NIR 
photometry) which is about 6 times larger than the one used in previous studies of dwarf galaxies. 
The NIR data are useful to study the L-Z relation because they are less affected by the ongoing 
star formation activities, and they also suffer significantly less from dust extinctions (both external 
and internal). Moreover, the Kg luminosity is a more reliable measurement of the stellar mass 
in these galaxies (Gil de Paz 2000; Perez-Gonzalez et al. 2003). Hence, the study of Lk^-Z 
relationship allows us to tie the observed variations in metal abundance to stellar mass, one of the 
most fundamental physical parameters of galaxies. However, there are only few works (Salzer et 
al. 2005; Mendes de Oliveira et al. 2006; Vaduvescu et al. 2007) studying the L-Z relation in the 
NIR due to the lack of photometric data. 

The paper is organized as following: Section 2 describes the sample of BCDs and data reduc- 
tions. In §3 we first derive the oxygen abundances for a subsample of BCD galaxies which have 
SDSS spectroscopical observations, and then we evaluate the B-, R- and Kg-hand L-Z diagrams. 
We compare our results with previous work, discuss the slopes of our L-Z relations and explore 
the possible origins of the scatter in the L-Z relation in §4. In the last section, a brief summary is 
present. 

2. Sample and Data Reduction 

There have been a number of criteria for selecting BCDs. These criteria are commonly based 
on the galaxy's morphological properties and its luminosity (Zwicky &; Zwicky 1971; Thuan & 
Martin 1981), although some definitions are based on their spectroscopic properties (Gallego et al. 
1996). BCDs classified by these different criteria may be confused with each other (Gil de Paz et al. 
2003). However, for our purpose, i.e. to determine the L-Z relationship, a clear and homogeneous 
sample of blue compact dwarf galaxies is necessary. 

Using a unified concept of BCD galaxy, Gil de Paz et al. (2003) compiled a BCD sam- 
ple (including 104 members) from several exploratory studies. The galaxies in their sample were 
selected by putting forward a new set of quantitative classification criteria. This new set of cri- 
teria is a combination of the galaxy's color, morphology and luminosity. Briefiy, a BCD galaxy 
has to fulfill the following three criteria (referring to blue, compact and dwarfs respectively): (1) 
Ai_B,pcak — A*i?,,pcak < 1, where fiB,pcak and /i_R,pcak are the peak surface brightness of B- and iZ-band, 
respectively, (2) /i_B,pcak < 22 mag arcsec"^, and (3) the absolute magnitude Mk > —21 mag [see 
Gil de Paz et al. (2003) for details]. Such a large and universal sample allows us to measure an 
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unbiased L-Z relation for BCDs. 

In this work, we select BCD galaxies from the Gil de Paz's sample, which is the largest well- 
defined sample of BCDs. Since not all galaxies have their metallicities measured, a galaxy is selected 
if (1) we can get its oxygen abundance from literature, or (2) its [Nil] A6584/Ha ratio was given in 
the original sample, or (3) we can obtain a spectrum from Sloan Digital Sky Survey (SDSS; Gunn 
et al. 1998; Blanton et al. 2003 ) Data Release Six (DR6). 

To derive the oxygen abundances for galaxies which have SDSS spectra, we need to measure 
the fluxes of emission-lines accurately. Thus we have to correct both the Galactic reddening and 
the underlying stellar absorption. The detailed method for the latter process can be found in 
Zhao et al. (2009; in preparation). We only give a brief description here. First, the observed 
spectra were dereddened for Galactic extinction, using the extinction coefficients from Schlegel et 
al. (1998) and the empirical extinction law from Cardelli et al. (1989). Then we model the stellar 
contribution through the stellar population synthesis code, STARLIGHT version 2.0 (Cid Fernandes 
et al. 2004). After subtracting the best-fit model spectrum from the observed one, we obtain the 
"pure" emission-line spectrum, from which we could measure accurate fluxes for all emission lines 
with the onedspec.splot task in IRAF0. The line flux errors are typically less than 5%. 



Finally, the observed line intensities were corrected for the intrinsic extinctions from the equa- 
tion 



/(H/5) F(H/5) 

where /(A) is the intrinsic line flux, F(X) is the observed line flux, c(H/3) is the logarithmic reddening 
parameter and /(A) is the redding function (Cardelli et al. 1989). For the intrinsic hydrogen line 
intensity ratios, we used the theoretical ratios from Brocklehurst (1971) at electron temperature and 
number density estimated from the observed [Oiii] A4959,5007/A4363 ratio and [Sii] A6717/A6731 
ratio, respectively (see §3.1 for detail). 

The photometric data for the optical bands are adopted from Gil de Paz et al. (2003). The NIR 
data used here are mainly taken from the Two Micron All Sky Survey (2MASS) extended source 
catalog (Jarrett et al. 2000). For several objects the Kg magnitudes are adopted from Noeske et 
al. (2003; 2005) and Vaduvescu et al. (2007). Because of the relatively shallower observations of 
2MASS, we prefer using the data from the later two sources to using the data from 2MASS when 
overlapping sources are found. 



^IRAF is distributed by the National Optical Astronomical Observatory, which is operated by the Association 
of Universities for Research in Astronomy (AURA), Inc., under a cooperative agreement with the National Science 
Foundation. 
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3. The L-Z Relationships 

3.1. The Metallicity 

Generally, the 0/H ratio is used to study the L-Z relationship since it is easily attained in 
the optical part of the spectrum. The preferred method for determining the oxygen abundance 
in galaxies is through electron temperature-sensitive lines (the so-called Tg method), such as the 
[O III] A4363 line (Aller 1984). The main problem is that normally the [O iii] A4363 line is very weak 
and only appears in very high excitation spectra. When the [O iii] A4363 line is absent, some other 
empirical methods, such as the ii23 method (Pagel et al. 1979), the P method (Pilyugin 2001) and 
the N2 method (Denicolo et al. 2002) have been employed frequently. These methods do not need 
any direct measurement of the electron temperature. 



3.1.1. Te Method 

We used the formula given by Izotov et al. (2006) to determine the abundances for galaxies 
which have good signal-to-noise ratios for the [O iii] A4363 line. Since the wavelength coverage of 
the SDSS spectra is 3800-9300 A, the [O ii] A3727 lines are not available for the low redshift (< 0.02) 
BCDs used here. Therefore, we used their equation (4), which makes use of the [On] A7320, 7331 
lines instead of the [On] A3727 line, to calculate the (O"*"/!!). The formulae is given the form of 

A7320 + A7331 2 487 

12 + log(0+/H+) = log — + 6.901 + 0.4831ogt - 0.013t + log(l - 3.48a;) (2) 

Hp t 

where t = 10~^Te[O iii], and x = 10~^net~*^'^. To calculate the (0^'^/H), we used the equation (5) 
in Izotov et al. (2006), as following 

12 + log(02+/H+) = i„g ^4959 + A5007 ^ ^ ^oo + 1:^ _ o.551ogt - O.OUt (3) 

Hp t 

The electron temperature (Tg) and number density (ng) were derived with the IRAF task nebu- 
lar. zones (De Robertis et al. 1987; Shaw & Dufour 1995) using the [Olii] A4363/(A4959 + A5007) 
ratio and the [S ii] A6717/A6731 ratio, respectively. This task is based on a five-level statistical 
equilibrium model and makes use of the latest collision strengths and radiative transition proba- 
bilities. For those we can not obtain the electron densities from the [S ii] line ratio, we assumed a 
density of 100 cm~^. This assumption almost does not affect our oxygen abundance determination, 
since the effect of temperature is much larger than that of electron density. For 12 out of 24 BCDs 
which have SDSS spectroscopical observations, we've measured the oxygen abundance using the Tg 
method, which is presented in Table 1. 

To study the L-Z relation, we also compiled oxygen abundance determined by the Tg method 
from the literature, as shown in Column (6) of Table 2. When more than one data point for a 
galaxy can be found, we use the error-weighted average value. For galaxies we can not obtain 
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errors of abundance from the original references, an equal weighted mean value is adopted. The 
final error is the bigger one of the two values: half of the difference of the range in values and the 
propagated error. 

3.1.2. N 2 Method 

The N2[=log([Nii] A6584/Ha)] empirical method is proposed by Denicolo et al. (2002), fol- 
lowing the earlier work by Storchi-Bergmann et al. (1994) and Raimann et al. (2000). They used a 
representative sample of Hll galaxies which have accurate oxygen abundances, plus photoionization 
models covering a wide range of abundances (7.2 < log(0/H) < 9.1), to calibrate the N2 estimator. 
As shown in Denicolo et al. (2002), the oxygen abundance and the N2 calibrator are well correlated 
(the correlation coefficient is 0.85), and the linear relation holds for the whole metallicity range, 
from the most metal-poor to the most metal-rich galaxies in their sample. Least-squares fits to the 
data simultaneously minimizing the errors in both axes, give 

12 -Mog(0/H) = 9.12-^0.73 xN2. (4) 

The oxygen abundance determined based on this method has a precision of ~ 0.2 dex. The 
N2 parameter has clear observational advantages for ranking metallicities in star-forming galaxies. 
Besides it being monotonic with log(0/H), it is also independent on reddening correction and flux 
calibration. What's more, for the BCDs having SDSS spectroscopic observations, the [O ii] A3727 
line, which is critical both for the i?23 method and the P method, is not available in the spectra. 
In addition to the above reasons, we have another reason for using N2 method other than the P 
or the i?23 method to derive the oxygen abundances for the BCD sample used in the current work. 
The [Nii]/Hq, ratio of about 47% galaxies were given by Gil de Paz et al. (2003), thus this method 
allows us to collect more data points than the other methods. The oxygen abundances determined 
based on the N2 method are listed in Table 1 (for the SDSS sample) and Column (7) of Table 2 
(for the entire sample). 

3.2. The L-Z Relations in the B- and i?-band 

In Table 2 we list the B and R magnitudes, along with their errors, which are taken from Gil 
de Paz et al. (2003). The foreground Galactic extinctions have been corrected for these galaxies 
using Ab values determined following Schlegel (1998) and the Galactic extinction law of Cardelli 
et al. (1989). The typical error in B magnitude is less than 0.05 mag, while it is around 0.1 mag 
in the i?-band. 

For the purpose of investigating the L-Z relation, we first display the broad band {B — R) color 
versus the metal abundance (derived by the Tg method) and the absolute B magnitude for our 
BCD sample, in the left and right panels of Figure 1, respectively. In the left panel, open circles 
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represent objects without a given measurement error of the abundance. One sees expected/famihar 
trends in both plots: the galaxies tend toward redder colors in progressing from lower to higher 
luminosity or metallicity, meanwhile with the presence of large scatters. These large scatters are 
mainly attributed to the internal redding and the measurement error. 

We present the Lb-Z and Lr-Z relations in Figures 2 and 3 respectively. Open symbols repre- 
sent the objects that the measurement errors of their abundances were not given in the referenced 
literature. However, the error will not affect our fitting results since we do not use it as a weighting. 

In the left panel of Figure 2 we show the 5-band absolute magnitude plotted versus metallicity, 
which was determined based on the Tg method. As we can see, a very good correlation is presented 
in the data. Only two galaxies, Mrk 328 and UM 533, which have the highest metallicities in the 
sample, might be discrepant from the trend. In order to check whether these two galaxies are really 
outliers, we calculated their deviations and found that both are less than 3(T. Thus, these two 
points don't really deviate from the trend. The apparent discrepancy may be caused by the fact 
that 1) the measurement error of the metallicity for UM 533 is as large as 0.29 dex, 2) BCDs with 
metallicity as high as > 8.6 are rare in our sample. A nonweighted least-squares linear fit, using a 
geometrical mean functional relationship (Isobe et al. 1990), to these 66 BCDs gives 

12 log(0/H)^^ = (4.53 ± 0.41) - (0.218 ± 0.038)Mij. (5) 

The fitted trend is shown as a solid line in Figure 2. The rms deviation of the data from the 
relationship is 0.25 dex in (0/H). The Spearman rank correlation coefficient (p, assessing how well 
an arbitrary monotonic function could describe the relationship between two variables) of the trend 
is -0.58 at a >4(T level of significance, which indicates an anti-correlation between Mb and O/H 
and means that the probability for the null hypothesis of no correlation between (0/H)t^ and Mb 
is ~ 0. The parameters of all fitting results presented in this section are summarized in Table 3. 

In order to illustrate the effect of using different methods to compute the abundance on the 
L-Z relation, we plot the N2-based L-Z relation in the right panel of the figures for each band 
considered. With a bit bigger scatter, a similar relationship also presents in the right panel of 
Figure 2. Excluding the labelled source. Pox 4 (its deviation is larger than 3cj), we obtain the 
following Lb-Z relation for the rest 73 BCDs from a geometrical mean fitting 

12 log(0/H)N2 = (4.18 ± 0.36) - (0.245 ± 0.037)Ms. (6) 

The rms deviation of the data from the relationship is 0.30 dex in (O/H). We can see that the 
Tg-based slope is in agreement with the N2-based slope within the uncertainty. However, as for the 
other two bands, the Tg-based slope is a bit smaller than the N2-based slope. We present a further 
discussion about this slope difference below. 

The situation changes a bit when coming to the longer wavelength photometric band, the R- 
band. The trend becomes shallower and the scatter is a bit smaller. While Pox 4 is still an outlier 
in this band. Nonweighted least-squares linear fits to these two data sets give 

12 log(0/H)y^ = (4.71 ± 0.36) - (0.199 ± 0.027)Mij. (7) 
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and 

12 + log(0/H)j^2 = (4-42 ± 0.32) - (0.221 ± 0.026)Mr. (8) 
for the Tg- and N2-based L^-Z relations respectively. 

3.3. The L-Z Relation in the K^-band 

The reason why we would like to arrive at an Lk^-Z relation is because such a relation should 
be more nearly representative of the more fundamental relationship between metal abundance and 
stellar mass. This is because the observed luminosity at the optical bands is dominated by recent 
star-formation activity rather than the stellar population that has accumulated over the galaxy's 
lifetime, while at the Kg-hand the light is sensitive to the bulk of the stellar content. 

The Li^^ relationship displayed in Figure 4 shows a change relative to the optical ones. The 
slope and scatter are much smaller than those of the Lb- and L^-Z relations, which can be also seen 
in Table 3. This is especially true for the Tg-based results: The rms scatter of the Lk^ -Z relation 
has about one third reduction comparing with those of the optical L-Z relations. We obtain the 
following relationships using nonweighted least-squares fits to the Tg-based and N2-based data sets 
respectively, 

12 + log(0/H)y^ = (4.97 ± 0.41) - (0.170 ± 0.024)Mii-^ . (9) 

and 

12 + log(0/H)N2 = (4.04 ± 0.38) - (0.221 ± 0.039)Mx^. (10) 

4. Discussion 

4.1. Comparison with the Literature 

While our study clearly demonstrates a good L-Z relation for BCD galaxies, the Mb — Z 
relation in Fig. 7 of HH99 is not obvious. Hopkins et al. (2002) also claimed that there is 
no correlation between oxygen abundance and total mass in their BCD sample. To explain the 
apparently contradictory results between these two studies and our work, we here compare our 
sample with HH99 and Hopinks' and discuss these results in the following. 

1) Sample size: The sample of HH99 only contains 13 BCDs (22 for Hopkins' sample) with ~4 
magnitude span, while our BCD sample is 3 — 5 times larger (see Table 3) and 2 magnitude wider. 
A small number statistics in sample size is difficult to avoid a selection bias. 

2) Uncertainties in determining 0/H: HH99 used a combination of oxygen abundances de- 
termined by (four) different methods for their BCD galaxies. Despite of a larger uncertainty in 
empirical-based 0/H comparing to Tg-based one, the inconsistency of the methods for determining 
0/H itself might cause an additional uncertainty in systematics. 
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3) In Hopkins et al. (2002), although all of their metallicities are derived with Tg-method, they 
used the width of H i emission-line profile at 20% of peak, W20, as an indicator of galaxy total mass. 
On the one hand, it requires radii and inclinations to calculate galaxy total masses from W2o- On 
the other hand, W20 also may not be a solid indicator for rotational velocity if turbulent motions 
are significant, as may be expected for BCD galaxies (e.g. Thuan & Martin 1981). 

We here used rather uniform data sets and systematically estimated 0/H in two different 
methods, and obtained the L-Z relationships in three bands {B, R &: Kg). Therefore, the results in 
the present work don't confiict with those in the literature, and BCDs do really show a well-defined 
L-Z relation. 

In Table 4 we list some relevant results of the L-Z relations compiled from the literature, 
and we plot them in Fig. 5. The slope of the Lb-Z relation for dwarf blue compact galaxies 
{Mb > —18; DBCGs; Shi et al. 2005) agrees with our results for BCDs within uncertainties. For 
dis {Mb > ~ —19), the slopes of the Tg-based Lb-Z relation are consistent with each other within 
the measurement errors (Skillman et al. 1989; Richer & McCall 1995; Lee et al. 2003; van Zee & 
Haynes 2006), and the mean value is —0.151, hence much shallower than ours (—0.218) for BCDs. 
However, the difference in slope between BCDs and dIs might not be truly significant, and could 
be mainly due to the small number statistics in sample sizes used in the study of the L-Z relation. 
The samples of dIs always include no more than 30 (i.e. 12 ~ 24) galaxies at the i?-band, and it is 
difficult to avoid a selection bias. It is possible that our results are less affected by sample selection 
given that the size of our sample is about 3 to 6 times larger. 

However, diversities do appear to be significant when we compare different slopes of the Lb- 
Z relation based on different empirical abundance calibrators. As shown in Table 4, all of the 
empirical calibrator-based slopes, except the absorption-free, KBG03-based one, for the sample 
{Mb > ~ —22) including massive galaxies in Salzer et al. (2005) are steeper than the Tg-based 
slopes for dIs, DBCGs and BCDs. But our N2-based slope is a bit steeper than or comparable with 
the absorption- free slopes from Salzer et al. (2005). At the same time, the different slopes derived 
with different abundance calibrators from Salzer et al. show large variations. These complicated 
situations might be the results of the fact that the empirical abundance calibrators are model- 
dependent, and thus have great complexity (see Kewley & Ellison 2008). This reminds us that we 
should be cautious when comparing slopes of the relation derived based on different abundance 
calibrators. 

Similar to the optical, the Tg-based slope of the NIR L-Z relation for BCDs is consistent with 
that for dIs (Mendes de Oliveira et al. 2006). While for the empirical abundance calibrator-based 
L-Z relation, our N2-based slope is in agreement with the (Te-|-i?23)-based slope given by Vaduvescu 
et al. (2007). We further discuss the various possible origins resulting in the differences in the slope 
and scatter in the L-Z relation below. 
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4.2. Slope of the L-Z relation 

4-2.1. The difference between T^- and N2-based slopes 

As shown in section 3.2 and 3.3, the slope derived with the N2-based metallicity is always some- 
what larger than that derived with the Tg-based metallicity for each photometric band (especially 
for the NIR //-.^relation), although they are consistent with each other within the uncertainty. In 
the following we discuss the possible origins of this difference. 

There are two effects which can cause the slope difference. One is the uncertainty in the 
abundance. The empirical abundance calibrator has larger measured uncertainty relative to the 
direct method, and thus might result in a slope difference. This effect is important when the sample 
is small. The other reason is the sample itself used to study the L-Z relation. In our work, the 
N2-based sample is a bit different from the Tg-based sample. These two effects could simultaneously 
affect the L-Z slope. 

To check this further, we use the same sample which both have direct and empirical metallicities 
to investigate the L-Z relations. In our sample, we find that there are 60 such galaxies in the optical 
and 37 such galaxies in the NIR. The slopes are -0.224 ± 0.038 and -0.223 ± 0.037, for the Tg- 
and N2-based Lb-Z relations respectively. These two slopes are very consistent. For the L^-Z 
relations, these two slopes are —0.204 ± 0.029 and —0.201 it 0.034, which are also in accordance 
with each other. This indicates that almost all of the slope change between the N2- and Tg-based 
optical L-Z relations is due to the sample issue. 

However, it is not the case for the NIR data. For the same sample, the N2-based slope is 
-0.204 ± 0.039, which is stiU larger than the Tg-based one, -0.171 ± 0.025. But the difference 
becomes smaller (see section 3.3). We find that the N2 method overestimates the metallicities for 
several objects at the high luminosity end when we plot these two data sets in the same figure. This 
can be also seen from Figure 6, in which the points with circles represent galaxies with M^g < —20 
mag. This overestimation (in fact, it results in a selection effect for a small sample), combining 
with the sample issue, causes the large difference between the T^- and N2-based slopes of the NIR 
relations. 

4.2.2. The NIR L-Z relation 

Now we return to the NIR data of BCDs. The slope of the L-Z relation for our BCD sample 
decreases monotonically as the wavelength of the photometric band used to construct it increases, 
which is in agreement with the result of Salzer et al. (2005). This trend can be seen directly 
through the green and red lines, for our and Salzer's results respectively, in Fig. 5. The dotted 
lines are the error-weighted linear fits to the data. Although our slopes have larger uncertainties 
(due to the small size of the sample), the tendency to decrease of the L-Z slope with wavelength is 
pronounced. 
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The variation of the slope of the L-Z' relation with wavelength might be expected simply due to 
stellar population issues. As shown by Bell & de Jong (2001), the characteristic colors of galaxies 
vary smoothly with M/L ratio, in the sense that higher M/L values are redder and lower M/L 
values are bluer. Color also correlates with luminosity, as shown in Fig. 1, with bluer galaxies 
being less luminous and redder ones more luminous. For our metal-poor systems, as an example, a 
galaxy with Mb = -12.0 and B-R = 0.0 will have a B - K color of ~2.0, and hence Mk « -14.0 
(see Table 4 of Bell & de Jong 2001). A galaxy with Mb = -18.0 and B - R = l.O will have 
B — K K, 3.0 and Mk ~ —21.0. Therefore, according to the models of Bell & de Jong (2001), one 
would expect a shallower slope in the iT-band L-Z relation, since in this example AM^ = 7.0 while 
AMs = 6.0 for the same two galaxies. 

Presumably, the effect of stellar population and the lower amounts of absorption will both be 
acting simultaneously to reduce the observed slope in the NIR L-Z relation. It is not clear which 
effect will be the dominant one (or whether they will be of roughly equal magnitude). To check this 
further, we try to investigate a more fundamental relation, MZ relation. We convert luminosity into 
stellar mass using the M/L ratio-color relation given by Bell & de Jong (2001), which is derived by 
using the Bruzual & Chariot (2003) model with a Salpter IMF (1955) and Z = OAZq. We adopt 
this sub-solar abundance model since our sample is metal-poor. The following three MZ relations 
are derived, 

12 + log(0/H)y_^ = (4.57 ± 0.34) - (0.408 ± 0.051)M^,ij, 

12 + log(0/H)y_^ = (4.57 ± 0.34) - (0.408 ± 0.051)M^,k, (11) 

12 + log(0/H)^^ = (4.72 ± 0.43) - (0.399 ± 0.054)M^,/^. 

where M^^b^ Af*,ij and M^^^k denote that the stellar mass is obtained by using the B-, R- and 
-fC-band luminosity respectively. From Eq. (11) we can see that the slopes of the MZ relations are 
very consistent. If we use the identical sample (i.e. same as the i^-band sample) for these three 
bands, the difference between these MZ slopes is even smaller, namely, 0.395 it 0.072 for B- and 
i?-band and 0.399 ± 0.054 for i^-band. These results suggest that the variation of M/L dominates 
the slope difference between the optical and NIR L-Z relations, while the absorption only has little 
effect, for our metal-poor galaxy sample. 

4.3. Scatter in the L-Z Relation 

The rms scatter of the Tg-based L-Z relation is much lower for the NIR data compared to the 
optical data. There are ~ 32% and ~ 26% reductions in the rms scatter of the Lks-Z relation 
relative to that of the Lb-Z and Lr-Z relation, respectively. This indicates that the NIR data 
should be more suitable for studying the L-Z relation. However, the rms scatter of the N2-based 
L-Z relation has little change (^ 10%). This is believed to be due to the relatively much larger 
uncertainties in the O/H derived with the empirical N2 method (see the following). 

Since the majority of the oxygen abundances are compiled from literature, our result may 
suffer from the uncertainty of aperture effect. We compare metallicities derived by different authors 
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(sources with more than one references for their metaUicities in Table 2). We find that the Tg-based 
oxygen abundances determined by different works agree with each other within 0.10 dex, except 
for a few points with large measurement errors. However, the N2-based abundances from different 
works appear to have larger scatter. Therefore, aperture effect may only have a minor contribution 
to the scatter in the Tg-based L-Z relation on condition that we've used an error-weighted mean 
value. While for the N2-based L-Z relation, this effect should be one main source causing the 
scatter. 

The N2-0/H relation (Eq. (4)), which itself has ^ 0.2 dex scatter, can also cause additional 
uncertainties in the N2-based L-Z relations. Generally, as shown in Fig. 6, the N2-based and Tg- 
based 0/H is uniformly distributed around the y = x line, but the standard deviation is as large as 
0.19 dex. For several galaxies, the deviation can be up to 0.4-0.6 dex. Therefore, the uncertainties 
introduced by the N2-method itself, along with the aperture effect, should be the major reason 
causing the deviant points in the N2-based L-Z plots. 

Given the fact that the typical errors are 0.05, 0.10 and 0.15 mag, in the -B-, R- and Kg- 
band photometry, respectively, the observed rms scatter of 1.15, 1.12 and 1.01 mag, in the 5-, R- 
and Ks-honA respectively, may be affected by some other parameters. Kobulnicky et al. (2003) 
searched for correlations between various parameters (e.g., color, size, EW of H/3, H/3 luminosity) 
and magnitude residuals (AM = the observed absolute magnitude — the derived magnitude from 
the relation using a given metallicity) in their derived Lb-Z relationship in order to investigate 
whether there are any "second-parameter effects" that might be responsible for increased scatter 
in the L-Z relation. They found no correlations for their local sample. 

Figure 7 shows the galaxy oxygen abundance, color, EWhq and Ha luminosity (from the left 
to the right) plotted versus AM from the best-fit linear Tg-based L-Z relation for BCDs. The 
bottom panels plot the full sample of 66 BCDs, while the upper two panels plot the subsample of 
39 galaxies which have Kg magnitude measured. Figures 7a and 7z indicate that there is a weak 
correlation between AM^ and metallicity. A Spearman rank correlation coefficient analysis gives 
p = 0.48 and p = 0.41 at a 2.9(t and 3.3(T level of significance, for Fig. 7a and Fig. 7i respectively. 
However, this correlation may be driven by the three points (in Fig. 7i) which have the highest and 
lowest metaUicities in this sample. To test this, we do the Spearman rank correlation coefficient 
analysis excluding these three points and obtain p = 0.32 at 2.5(T confidence. This suggests the 
scatter of Lb-Z relation does really correlate with the metallicity, although the correlation is weak. 

Meanwhile, both of Figures Id and 7/ show anti-correlations between AM^ and Lho, indicating 
that this correlation really exists. A Spearman rank correlation coefficient analysis gives p = 0.55 at 
a > 4(7 level of significance for Fig. 7/. We got a similar result (p = 0.53 at > Aa confidence) when 
excluding the two points with the largest deviations (one has a huge metallicity uncertainty and the 
other doesn't have a metallicity uncertainty) and the point with the largest Lho- Therefore, these 
two correlations are not driven by some suspect points but reasonable in the sense that (1) a higher 
metallicity might possibly cause a larger internal extinction; (2) higher Ha luminosity means higher 
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blue luminosity. But the large scatters in these two correlations indicate that the observed scatter 
may represent some fundamental cosmic scatter (e.g. time of formation of galaxies, age dispersion, 
etc; see Calura et al. 2009) in the L-Z relation. However, no obvious correlation between the {B-R) 
or EWho and AM can be seen in Figures 76, 7/, 7j, 7g and 7k. By comparing Figure 7c with Figure 
7k, one can find that the correlation appearing in Figure 7c is due to the small number of statistics 
involved in sample. 

For the ifg-band, there only exists one possible correlation between AM;<-^and metallicity, and 
we can not draw a firm conclusion because of the lack of data at the low metallicity end. To check 
this, we compare the same data set of the 5-band with the ifg-band, and find a similar trend for 
all of three samples (see Figures 7a, 7e and 7i). Therefore, the correlation between the L-Z relation 
residuals and metallicity is true, and is not due to the sample issues. However, a Spearman rank 
correlation coefficient analysis indicate that the correlation between AM^^^and Z \s much weaker 
than that in 5-band. The correlation coefficient is 0.34 and only at a 2.1(7 level of significance. 
This is because that the NIR data suffer much less extinctions than the optical data. In contrast 
to the i?-band, the seeming correlation between AM/^^ and Lho, as shown in Figure 7/i, is caused 
by the sample issue. This result is very reasonable since the current star forming activities should 
have much less effect on Lk,, than on Lb- 



5. Summary 

By cross-correlating the blue compact dwarf galaxy sample of Gil de Paz et al. (2003) with 
the SDSS DR6, we derived and compiled oxygen abundances for a subsample containing ~ 70 BCD 
galaxies using both the direct (Tg) and empirical (N2) methods. The uniformity of the sample, 
combining with its relatively large size, makes this an excellent data set to systematically study 
the L-Z relations for BCDs for the first time. We investigate the L-Z relations in three different 
photometric bands {B, R, and Kg) and explore the effects of using different abundance calibrators 
on the resulting L-Z relations. The main results are summarized bellow: 

(1) We present L-Z relations for all three bands, from optical to near infrared, for these BCD 
galaxies. Our study reconciles apparently contradictory results obtained by different authors. We 
find that the slope of L-Z relation calibrated with the empirical N2-method, for each photometric 
band, is consistent with that evaluated using the Tg-based metallicities within the uncertainty. 

(2) We confirm the similar result for BCDs here that the slope of the L-Z relation becomes 
shallower with longer wavelength of the photometric band used, as has been found in larger samples 
that include more luminous systems (Salzer et al. 2005). 

(3) The correlation between the L-Z relation residuals and Hq luminosities in the i?-band, 
combined with the fact that no such correlation exists in the i^'^-band, indicates that the starburst 
activities might be a cause of the large scatter in the Lb-Z relation. Meanwhile, the weak correlation 
between the L-Z relation residuals and metallicity might suggest that the internal-absorption is 
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another possible factor contributing to the scatter in the L-Z relationship. 
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Table 1. Electron temperatures, number densities and oxygen abundances of BCDs in SDSS 
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12+log{0/H) 


Object Name 


(K) 


(cm~'^) 
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HS 0822+3542 
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7.41 


7.49 


7.42*^ 


HS 1400+3927 
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8.11 


8.06 




HS 1440+4302 


12544. 


49. 


8.10 


8.05 


8.07= 


HS 1609+4827 








8.37 




Haro 3 
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57. 


8.56 


8.46 


8.29*' , 8.37= , 8.37'* 


II Zw 71 








8.40 




Mrk 1313 


11921. 


26. 


8.17 


8.08 


8.19*= , 8.34= 


Mrk 1416 
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7.82 
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Mrk 1418 








8.48 




Mrk 1423 








8.52 




Mrk 1480 








8.03 




Mrk 1481 








8.11 




Mrk 178 








7.86 




Mrk 409 








8.81 




BBS 1054+504 


10610. 


100.=^ 


8.26 


8.34 




BBS 1147+520 








8.19 




BBS 1428+457 


9697. 


95. 


8.44 


8.18 


8.40*^ 


BBS 1533+574 


11690. 


21. 


8.14 


8.13 


8.00= 


UGCA 184 .. . 


12825. 


100.=" 


8.04 


7.86 


8.04*= , 8.00*" 


UM 439 


12099. 


24. 


8.09 


8.06 
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100.^ 


8.27 


8.40 




UM 456A .... 








8.28 
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8.31 




VCC 0130 .... 








8.40 





''ne=100. cm ^ is an assumed nominal value. 



References. — ^ Izotov et al. (2006); = Nagao et al. (2006); <* Steel et al. 
(1996); = Masegosa et al. (1994); * Kobulnicky & Skillman (1996). 
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Table 2. Sample and properties of BCDs used to study the L-Z relationships 
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29 


Mrk 0005 


30 


60 


15 


13 


± 





04 


14 


56 


± 





13 


12.68 


8.06 


+ 


04 


8 


19 


3 


5 


29, 


34 


Mrk 108 


31 


69 


15 


15 


± 





03 


14 


66 


± 





08 


13.76 + 0.20 


7.96 


+ 


02 


8 


03 


1 


15, 16 




29 


Mrk 1313 


32 


51 


16 


02 


± 





03 


15 


50 


± 





06 


13.63 + 0.27 


8.22 


+ 


10 


8 


05 


1 


6, 7, 18 


6, 7, 


29 


Mrk 1329 


31 


02 


14 


08 







03 


13 


38 







09 


11.93 


8.26 


+ 


02 


7 


95 


2 


5, 13 


28, 


29 


Mrk 1416 


32 


64 


16 


32 


± 





03 


15 


75 


± 





02 




7.84 


+ 


02 


7 


94 




5, 6 


6, 


32 


Mrk 1418 


30 


28 


13 


86 


± 





03 


12 


61 


± 





19 


11.604 + 0.061 








8 


48 


1 






6 


Mrk 1423 


31 


55 


14 


90 







03 


13 


52 







10 


12.65 + 0.11 








8 


52 


1 






6 


Mrk 1434 


32 


62 


16 


77 


± 





04 


16 


43 


± 





13 




7.79 


+ 


01 


7 


69 




5 




29 


Mrk 1450 


30 


83 


15 


75 


± 





05 


15 


09 


± 





09 




7.96 


+ 


02 


7 


93 




5 




35 


Mrk 1480 


32 


18 


16 


17 


± 





03 


15 


56 


± 





08 


13.66 + 0.18 


8.04 


+ 


05 


8 


03 


1 


7 


6 


, 7 


Mrk 1481 


32 


18 


16 


19 


± 





03 


15 


57 


± 





06 










8 


11 








6 


Mrk 178 


28 


11 


14 


15 


± 





04 


13 


60 


± 





09 


11.74 


7.92 


+ 


02 


7 


88 


4 


5, 12 


6, 


29 


Mrk 209 


28 


82 


14 


15 


± 





03 


13 


94 


± 





13 


12.70 


7.76 


+ 


01 


7 


66 


4 


5 




29 


Mrk 324 


32 


01 


15 


17 


± 





03 


14 


60 


± 





13 


12.50 + 0.06 








8 


18 


3 






29 


Mrk 328 


31 


75 


14 


93 


± 





05 


14 


18 


± 





09 


12.10 + 0.08 


8 


.66 




8 


64 


1 


8 




29 


Mrk 409 


31 


64 


14 


37 


± 





03 


13 


33 


± 





05 


11.17 + 0.03 








8 


80 


1 






6 


Mrk 450 


30 


57 


14 


44 


± 





05 


13 


65 


± 





10 


11.97 


8.12 


+ 


02 


8 


09 


3 


5 


28, 


29 


Mrk 475 


29 


93 


16 


20 


± 





03 


15 


67 


± 





07 




7.93 


+ 


02 


7 


89 




5 




35 


Mrk 600 


30 


82 


14 


85 


± 





03 


14 


82 


± 





15 


13.11 


7.94 


+ 


06 


7 


84 


3 


5, 12, 17 


29, 


34 


Mrk 67 


30 


77 


16 


10 


± 





03 


15 


34 


± 





08 


14.38 + 0.26 


8.08 


+ 


08 


7 


90 


1 


5, 11 


28, 


29 


Mrk 709 


30 


98 


16 


32 


± 





03 


15 


65 


± 





02 


14.14 + 0.27 


7.68 


+ 


04 


8 


38 


1 


18 




29 


Mrk 86 


29 


20 


12 


07 


± 





03 


11 


49 


± 





11 


9.13 + 0.03 








8 


53 


1 






29 


Mrk 900 


31 


37 


14 


17 


± 





03 


13 


56 


± 





12 


11.38 + 0.10 


8.07 


+ 


03 


8 


48 


1 


19 




29 


Mrk 996 


31 


88 


15 


01 


± 





03 


14 


08 


± 





14 


11.85 








8 


21 


3 






29 


NGC 1522.... 


30 


13 


14 


03 


± 





01 


13 


22 


± 





06 


11.38 + 0.10 


8.07 


+ 


05 






1 


18 






NGC 1705.... 


28 


54 


13 


09 


± 





01 


12 


19 


± 





06 


10.52 + 0.06 


8.21 


+ 


05 






1 


20 






NGC 2915.... 


27 


78 


11 


93 


± 





01 


10 


96 


± 





06 


9.83 + 0.06 








8 


40 


1 






36 
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Table 2 — Continued 



Object Name 
(1) 


DM^' 
(mag) 
(2) 




(mag 
(3) 


) 






R" 
(mag 
(4) 


) 






(mag) 
(5) 




12+log(0/H) 
Te N2 
(6) (7) 


(8) 


Reference 
(9) 


(10) 


NGC 3125 


29 


84 


13 


05 


± 





01 


12 


25 







06 


10 


52 + 0.05 


8 


29 


+ 


10 


8 


01 


1 


9, 12, 18, 21 




29 


NGC 4861 


30 


50 


12 


68 


± 





03 


11 


91 


± 





09 


11 


77 + 0.11 


8 


00 


+ 


01 


7 


90 


1 


5, 12, 22 


29, 


32 


Pox 186 


30 


93 


17 


73 


± 





03 


17 


39 


± 





01 






7 


74 


+ 


01 








12, 23 






Pox 4 


33 


45 


15 


27 


± 





04 


14 


88 


± 





10 




13.52 


7 


98 


+ 


06 


7 


67 


2 


12, 21, 23 




29 


SBS 0940+544C 


31 


91 


17 


18 


± 





04 


17 


05 


± 





10 






7 


43 


+ 


02 


7 


50 




5 




29 


SBS 1054+504 . . 


31 


52 


16 


08 


± 





04 


15 


46 


± 





15 


13 


13 + 0.14 




8 


.26 




8 


34 


1 


6 




6 


SBS 1147+520. . 


31 


39 


16 


95 


± 





05 


15 


98 


± 





09 














8 


19 








6 


SBS 1331+493 . . 


29 


98 


14 


87 


± 





03 


14 


16 


± 





17 






7 


78 


+ 


02 


7 


89 




5 


29, 


35 


SBS 1415+437 . . 


30 


02 


15 


43 


± 





03 


14 


77 


± 





12 






7 


59 


+ 


01 


7 


76 




5, 24 


29, 


34 


SBS 1428+457 . . 


32 


74 


15 


42 


± 





05 


14 


67 


± 





09 


13 


15 + 0.14 


8 


42 


+ 


05 


8 


18 


1 


6, 7 




6 


SBS 1533+574. . 


33 


48 


16 


02 


± 





10 


15 


30 







10 


13 


74 + 0.21 


8 


07 


+ 


07 


8 


11 


1 


5, 6 


6, 


32 


Tol 002 


29 


25 


14 


06 


± 





01 


13 


34 


± 





06 


12 


47 + 0.13 


7 


98 


+ 


02 


8 


15 


1 


12, 18 




29 


Tol 1345-420 .... 


32 


54 


15 


87 


± 





03 


15 


00 


± 





03 


14 


11 + 0.22 


7 


99 


+ 


02 


7 


85 


1 


12, 18 




29 


Tol 1434+032 . . . 


31 


66 


15 


86 


± 





03 


15 


42 


± 





06 






7 


88 


+ 


04 


7 


87 




7 




7 


Tol 17 


32 


19 


15 


99 


± 





05 


15 


05 


± 





06 






7 


94 


+ 


05 


8 


17 




18 




29 


Tol 35 


32 


16 


14 


17 


± 





03 


13 


22 


± 





05 


11 


73 + 0.08 


8 


19 


+ 


01 


8 


08 


1 


12, 18 




29 


Tol 65 


32 


78 


17 


26 


± 





04 


16 


84 


± 





03 




15.99 


7 


50 


+ 


01 


7 


51 


2 


5, 12, 18 


29, 


37 


Tol 85 


33 


47 


16 


51 


± 





03 


15 


99 


± 





01 






8 


04 


+ 


01 








18 






UCM 1612+1308 


33 


63 


17 


21 


+ 





14 


17 


31 


+ 





07 






8 


17 


+ 


03 


7 


96 




25 




29 


UGCA 184 


31 


81 


15 


99 


± 





04 


15 


80 


± 





10 


14 


17 + 0.22 


8 


03 


+ 


01 


7 


89 


1 


5, 6, 7 


6, 


35 


UGCA 412 


33 


07 


15 


55 


± 





05 


14 


62 


± 





09 








8 


.12 




8 


38 




8 




8 


UM 133 


31 


91 


15 


41 


± 





03 


14 


51 


± 





14 






7 


69 


+ 


02 


7 


78 




5 




28 


UM 323 


32 


26 


16 


09 


± 





04 


15 


24 


± 





08 






7 


96 


+ 


04 


8 


00 




7 




7 


UM 382 


33 


53 


18 


20 


± 





04 


17 


83 


± 





22 






7 


82 


+ 


03 








26 






UM 408 


33 


58 


17 


46 


± 





16 


16 


70 


± 





12 






7 


74 


+ 


05 


7 


90 




18, 27 




27 


UM 439 


30 


73 


14 


77 


± 





03 


14 


09 


± 





06 






8 


08 


+ 


03 


7 


93 




5, 6, 12 


6, 


28 


UM 452 


31 


45 


15 


25 


± 





03 


14 


07 


± 





11 


12 


92 + 0.17 




8 


.27 




8 


40 


1 


6 


6, 


29 


UM 455 


33 


64 


17 


02 


± 





03 


16 


26 


± 





01 






7 


74 


+ 


02 








18 






UM 456A 


31 


95 


16 


71 


± 





03 


16 


05 


± 





01 














8 


21 








6 


UM 491 


32 


20 


15 


54 


± 





03 


14 


95 


± 





07 


13 


43 + 0.21 










8 


31 


1 






6 


UM 533 


30 


34 


14 


63 


± 





04 


13 


64 


± 





11 






8 


55 


+ 


29 


8 


32 




7 




7 


VCC 0130 


31 


02 


17 


05 


± 





05 


16 


27 


± 





01 














8 


28 








6 


VCC 0459 


31 


02 


14 


95 


± 





05 


14 


13 


± 





01 




12.45 


8 


27 


+ 


09 


8 


28 


4 


4 




4 


VCC 0655 


31 


02 


13 


32 


± 





04 


12 


12 


± 





03 


10 


48 + 0.02 










8 


72 


1 






38 


VCC 0848 


31 


02 


15 


03 


± 





05 


14 


10 


± 





01 




12.71 


8 


06 


+ 


12 


8 


27 


4 


5 




38 


VII Zw 403 


28 


41 


14 


11 


± 





04 


13 


58 


± 





11 




12.48 


7 


70 


+ 


01 


7 


84 


4 


5 




32 



Note. — Col. (1): Galaxy name. Col. (2): Distance moduli. Col. (3) - (5): B, R and Ks-hand apparent magnitudes, 
respectively. Col. (6): Oxygen abundances determined by the Te method. Col (7): Oxygen abundances determined by the N2 
method. Col. (8): References for the Kg magnitudes. Col. (9): References for the metal abundances derived by the Te method. 
Col. (10): References for the N2 ratios. 

''AH data are from Gil de Paz ct al. (2003), and all of the magnitudes have been corrected for the Galactic extinctions. 

''The Galactic extinctions are not corrected for the magnitudes from 2MASS, while the corrections have been applied for the 
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magnitudes from the other three sources. 

References. — (1) Jarrett et al. 2000; (2) Noeske et al. 2003; (3) Noeske et al. 2005; (4) Vaduvescu et al. 2007; (5) Nagao et al. 
2006; (6) This work; (7) Izotov et al. 2006; (8) Shi et al. 2005; (9) Mas-Hesse & Kunth 1999; (10) Steel et al. 1996; (11) Garnett 
1990; (12) Kobulnicky & Skillman 1996; (13) Guseva et al. 2000; (14) French 1980; (15) AUoin et al. 1979; (16) van Zee et al. 
1998; (17) Augarde et al. 1990; (18) Masegosa et al. 1994; (19) van Zee & Haynes 2006; (20) Lee & Skillman 2004; (21) Vacca 
& Conti 1992; (22) Koubulnicky & Skillman 1998; (23) Guseva et al. 2007; (24) Thuan et al. 1999; (25) Rego et al. 1998; (26) 
Kniazev et al. 2001; (27) Pustilnik et al. 2002; (28) Izotov & Thuan 2004; (29) Gil de Paz et al. 2003; (30) Melbourne et al. 2004; 
(31) Guseva et al. 2003; (32) Izotov et al. 1997; (33) Jansen et al. 2000; (34) Izotov et al. 1998; (35) Izotov et al. 1994; (36) Lee 
et al. 2003; (37) Izotov et al. 2004; (38) Vi'chez & Iglesias-Paramo 2003. 
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Table 3. Results of linear fits to the L-Z relations 



Method 


Band 


No. of Galaxies 


Intercept 


Slope 


rms 


P 


Significance 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 




B 


66 


4.53 ±0.41 


-0.218 ±0.038 


0.25 


-0.58 


>4o- 




R 


66 


4.71 ±0.36 


-0.199 ±0.027 


0.23 


-0.63 


>4(T 




Ks 


39 


4.97 ±0.41 


-0.170 ±0.024 


0.17 


-0.68 


>4cr 




B 


73 


4.18 ±0.36 


-0.245 ± 0.037 


0.30 


-0.53 


>4o- 


N2 


R 


73 


4.42 ± 0.32 


-0.221 ±0.026 


0.29 


-0.60 


>4(T 




Ks 


47 


4.04 ± 0.38 


-0.221 ±0.039 


0.27 


-0.56 


3.7o- 



Note. — Col. (1): Method for determining the oxygen abundance. Col (2): Wavelength 
band. Col (3): Number of galaxies used in the fit. Col (4): L-Z relation intercept. Col (5): L-Z 
relation slope. Col. (6): rms error in log(0/H). Col. (7): Spearman rank correlation coefficient. 
Col (8): Level of significance for the Spearman rank correlation coefficient. 
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Table 4. Luminosity-metallicity relations for star-forming galaxies 



Sample 




Method 


Intercept 


Slope 


rms 


Reference 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


iJ-band 


19 dis 




Te 


5.50 


—0.153 


0.16 


SKH89 


12 dis 




Te 


5.67 ±0.48 


-0.147 ±0.029 


0.09 


RM95 


22 dIs 




Te 


5.59 ±0.54 


-0.153 ±0.025 


0.18 


LMK03 


24 BCGs Mb > 


-18 


Te 


5.50 ±0.81 


-0.16 ±0.05 


> 0.2 


SKL05 


21 dis 




Te 


5.65 ±0.17 


-0.149 ±0.011 


0.15 


ZH06 


66 BCDs 




Te 


4.53 ±0.41 


-0.218 ±0.038 


0.25 


Present work 


29 dis 




Te+P 


5.80 ±0.17 


-0.139 ±0.011 


~ 0.15 


PVC04 


24 BCGs Mb > 


-18 


N2 


5.16 ±0.63 


-0.19 ±0.05 


> 0.2 


SKL05 


73 BCDs 




N2 


4.18 ±0.36 


-0.245 ± 0.037 


0.30 


Present work 


54 Spirals 




P 


6.93 ±0.37 


-0.079 ±0.018 


~ 0.2 


PVC04 


~ 53000 SFGs 




SDSS R23 


5.238 ±0.018 


-0.185 ±0.001 


> 0.3 


THK04 


48 BCGs Mb < 


-18 


Te 


7.37 ±0.53 


-0.05 ±0.03 


> 0.2 


SKL05 


48 BCGs Mb < 


-18 


N2 


7.06 ± 0.63 


-0.08 ±0.03 


> 0.2 


SKL05 


765 ELGs 




EP84 R23 


3.75 ±0.05 


-0.280 ±0.003 


0.29 


SLM05 


765 ELGs 




KBG03 ^23 


4.56 ±0.05 


-0.222 ± 0.003 


0.25 


SLM05 


765 ELGs 




SDSS ^23 


3.96 ±0.05 


-0.271 ±0.003 


0.28 


SLM05 


765 ELGs 




EP84 _R23 


3.75 ±0.05 


-0.253 ±0.002 


0.26 


SLM05'' 


765 ELGs 




KBG03 R23 


4.56 ±0.05 


-0.200 ± 0.002 


0.22 


SLM05'' 


765 ELGs 




SDSS ^23 


3.96 ±0.05 


-0.247 ±0.002 


0.25 


SLM05=' 


R-band 


765 ELGs 




EP84 _R23 


3.42 ±0.05 


-0.266 ±0.002 


0.27 


SLM05'" 


765 ELGs 




KBG03 R23 


4.28 ±0.06 


-0.211 ±0.003 


0.23 


SLM05'' 


765 ELGs 




SDSS R23 


3.64 ±0.05 


-0.259 ±0.002 


0.26 


SLM05'' 


66 BCDs 




Te 


4.71 ±0.36 


-0.199 ±0.027 


0.23 


Present work 


73 BCDs 




N2 


4.42 ±0.32 


-0.221 ±0.026 


0.29 


Present work 


/^3-band 


370 ELGs 




EP84 _R23 


3.92 ±0.09 


-0.212 ±0.003 


0.24 


SLM05 


370 ELGs 




KBG03 ^23 


4.03 ±0.10 


-0.195 ±0.004 


0.23 


SLM05 


370 ELGs 




SDSS ^23 


4.85 ±0.09 


-0.173 ±0.003 


0.18 


SLM05 


27 dis 




Te 


5.37 ±0.20 


-0.151 ±0.013 


0.14 


OTC06== 


25 dis 




Te + R23 


5.58 ±0.24 


-0.141 ±0.015 


0.10 


VMR07 


14 BCDs 




Te + R23 


4.21 ±0.54 


-0.224 ±0.030 


0.11 


VMR07<i 


39 BCDs 




Te 


4.97 ±0.41 


-0.170 ±0.024 


0.17 


Present work 


47 BCDs 




N2 


4.04 ± 0.38 


-0.221 ±0.039 


0.27 


Present work 



Note. — Col. (1): Sample: dis, dwarf irregular galaxies; SFGs, star-forming galaxies; BCGs, 
blue compact galaxies; ELGs, emission line galaxies. Col. (2): Method for determining the oxygen 
abundance. Col (3): L-Z relation intercept. Col (4): L-Z relation slope. Col. (5): rms error in 
log(0/H). Col. (6): Reference. 

^These results are corrected for the internal absorption in luminosity. 

'^Results of the L-Z relation for V-band. 

'^Results re-estimated by us using a geometrical mean fit to their 27 data points (excluding the two 
points with metallicity derived by the R23 method). 

"^It is worth noting that the definition of BCDs used in their work is different from that in the 
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present work. 

References. — (SKH89) Skillman et al. 1989; (RM95) Richer & McCall 1995; (LMK03) Lee et 
al. 2003; (SKL05) Shi et al. 2005; (ZH06) van Zee & Haynes 2006; (PVC04) Pilyugin et al. 2004; 
(THK04) Tremonti et al. 2004; {SLM05) Salzcr et al. 2005; (OTC06) Mendes de Olivcira et al. 2006; 
(VMR07) Vaduvescu et al. 2007. 
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7.5 8.0 8.5 -18 -17 -16 -15 -14 -13 -12 

12 + log(0/H)T^ Mb 



Fig. 1. — Left panel: The (B — R) color vs. metallicity diagram. Open circles show objects with 
errors of the abundances unavailable. Right panel: The {B — R) color vs. the 5-band absolute mag- 
nitude of the BCD galaxies. The familiar trends (i.e. objects with higher luminosity or metallicity 
having redder color) are clearly present in these two plots. 
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Fig. 2. — 5-band L-Z relationships for the BCDs. The objects without measurement errors in 
abundances and/or magnitudes are plotted with open symbols. Larger symbols show objects which 
are not used to construct the L-Zrelations. The solid lines represent geometrical mean, least-squares 
linear fits to the data. Lejt panel: Tg-based results. Right panel: N2-based results. 
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Fig. 3. — Same as Fig. 2, but for the i?-band. 
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Fig. 4. — Same as Fig. 2, but for the Ks-hand 
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-0.25 



SKHB9, RM95, SKL05, ZH06, OTC06 
Present work 
Present work 

□ PVC04, SKL05, VMR07 
■ SKL05 

PVC04, THK04, SKL05 

□ PVC04, SKL05, VMR07 




> SLM05: EP84 

< SLM05: KBG03 

E SLM05: SDSS 

Filled: Te method 

Open: empirical method 

Green: dwarf galaxies 

Red: giant/mixed galaxy sample 



3.6 



3.8 



4.0 

log A (A) 



4.2 



4.4 



Fig. 5. — Compiled slopes of Z-Z relations. For clarity, some points for B- and Kg-hand are shifted 
a bit along the x-axis. The dotted lines are the results of error-weighted linear fits to the data. The 
names of the references are the same as those in Table 4. 
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7.5 8.0 8.5 

12 + log(0/H)T 



Fig. 6. — Comparison of Tg-based (abscissa) and N2-based (ordinate) 0/H. Points with circles 
represent galaxies with Mk^ ^ —20.0 mag. The solid line is a reference line for the case when the 
two quantities are the same. 
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Fig. 7. — Galaxy oxygen abundance, color, EWh^ and Ha luminosity (from the left to the right) 
vs. luminosity residuals, AM, from the best-fit linear Tg-based L-Z relation. The bottom panels 
plot the full sample of a total of 66 BCDs, while the upper two panels plot the subsample of 39 
galaxies which have magnitudes. Panels (a), (e) and {i) indicate that there is a weak correlation 
between AM and metallicity. While no correlation between AM and [B — R) color or EWho can 
be investigated [see panels (&), (/), (j), (5) and {k)]. By comparing panel (c) with (A;), one can find 
that the correlation present in panel (c) is due to selection bias. This bias is also present in panel 
[K). Meanwhile, both of panels {d) and (Z) show correlations between AM^ and Lhq, indicating 
that this correlation really exists. 



